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CALCULATION  xir  GOLDSISEf  FACTORS 


In  his  formulation  of  the  vortex  theory  of  screw  propellers, 

S.  Goldstein^)* derives  the  following  expression  for  the  distribution 
of  circulation  about  each  blade  section: 
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where  the  coefficients  satisfy  the  following  infinite  system  of 
linear  equations: 
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The  Goldstein  factor  K  is  defined  in  terms  of  the  circulation 
by  the  relation 

K  = 

*  References  axe  listed  on  page  17* 


In  the  preceding  equations,  p  denotes  the  number  of  blades,  T*  is 
the  circulation,  id  is  the  angular  velocity  and  v  is  the  velocity  of  advance 
at  the  propeller,  v  ic  the  velocity  of  the  screv  surface  in  the  direction 
of  its  axis,  y-  represents  vr/u  vhere  r  is  the  distance  from  the  axis 
of  rotation,  and  fL  is  the  value  of  when  r  ■  5,  the  radius  of  the 
propeller.  The  symbol  ■$>*.+£)  ^  ^  represents 

Golds  tain’s  function,  which  is  defined- by  the  relation 

T  to  =  5  (**•)  (4) 

i,>0  '  tj-rC  *>v 

vhere  8^  n(ix)  is  a  Xcomel  function^2),  and  K^(x)  is  the  modified  Bessel 
function  of  the  second  kind,  of  order  n.  Similarly,  I^x)  represents  the 
modified  Bessel  function  of  the  first  kind,  of  order  n.  Primes  are  affixed 
to  these  functions  to  denote  differentiation  with  respect  to  the  argument. 

For  large  values  of  h>  Goldstein  showed  that  the  individual 
terms  of  the  equations  comprising  the  infinite  system  (2)  are  closely 
approximated  by  the  corresponding  terms  of  the  infinite  system 
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and  he  succeeded  in  solving  this  system  exactly,  obtaining 


a*  - - zt- 


(A»v)/ 


Wv 


In  the  solution  of  AML  Problem  42  -  54,  equation  (l)  as  modified 
by  the  substitution  of  a^*  for  was  evaluated  on  UNI7AC  for  values 
of  ranging  from  0.25  to  6;  p  »  3,  4,  5>  6;  ana  extending 


(6) 


2 


from  0.2  to  0*975*  Subsequent  desk-machine  calculations  have-  revealed 
that  the  use  of  V  in  place  of  3eads  to  errors  in  the  calculated, 
circulation  that  are  unacceptably  large  when  fa^  is  less  than  2. 

Accordingly,  further  research  was  requested,  requiring  an  investi¬ 
gation  of  procedures  for  solving  the  infinite  system  (2)  for  a  sufficient 
number  of  the  successive  unknowns  aQ  (m  =  0,  1,  2,  to  an  adequate 
degree  of  accuracy  to  yield  reliable  values  of  the  circulation. 

The  purpose  of  this  report  is"  to  set  forth  in  detail  an  effective 
procedure  for  solving  the  infinite  system  (2)  for  the  leading  a^'s 
when  fa  exceeds  unity.  This  limitation  of  the  size  of  fa  is 
necessary  to  insure  a  satisfactory  rate  of  convergence  of  the  successive 
approximations  to  the  unknowns  a^.  Fortunately,  the  practically  important 
range  of  fao  coincides  with  the  one  that  we  are  considering. 

For  the  sake  of  clarity,  we  shall  illustrate  the  proposed  technique 
by  applying  it  to  the  calculation  of  the  leading  unknowns  aQ,  a^,  .  . 
corresponding  to  the  particular  values  =1*5  and  p  =  3* 

The  coefficients  occuring  in  the  equations  of  system  (2)  are 
evaluated  by  a  variety  of  methods.  When  the  arguments  and  orders  of 
the  Bessel  functions  IQ  and  and  their  derivatives  3^'  and  K^'  are  less 
than  20  the  functional  values  can  be  obtained  from  existing  tables(3>^). 
Outside  this  range  of  arguments  and  orders,  the  following  formulas ( 5) 
are  effective. 
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Goldstein  function  T1>n(x)  is  calculable  by  the  series 


T->;  =  '  ifcr  ~  wkf^i  -  '■■») 


•/n. 


provided  n  is  not  an  even  integer. 

When  n  is  an  even  integer, 

T  (*■)  =  I  -  *£  +  ^  +  ■  • 

I,*  X  *  4 
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An  alternative  series  for  ^  is 

T  K“J  ~  +  vJrJ  +  .  tvfoL  +  •  • 


TV 
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vhere 
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The  values  of  Tj  (x)  cen  "be  calculated  by  means  of  6eri.es  obtained 
l,n 

by  term-by-term  differentiation  of  the  preceding  series  for  T.  (x). 

l^n 

With  the  aid  of  such  formulas,  the  leading  coefficients  of  the 
system  (2)  vere  computed  when  yLC  -  1.5  and  p  =  3.  Jhe  resulting  system 
is  shown  on  the  following  page. 

The  coefficient  of  in  the  nth  equation  of  this  infinite  system 
is  closely  approximated  by 


T  o±j£/x 


and  the  constant  term  is  approximated  by 
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Indeed,  these  approximations  to  the  coefficients  and  constant  terms 
correspond  to  the  simplified  system  (5). 

Our  procedure  for  solving  the  system  (12)  is  as  follows.  We 
initially  restrict  our  attention  to  the  first  equation  of  (12)  and 
assume  that 

a.  =  a.*  = _ ai. _ _ 

for  all  m«S£  1.  From  equation  (5)  we  infer  that 
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(14) 


(15) 


p  yyt.-o 

and  consequently  the  first  equation  of  system  (12)  can  be  replaced  by 
the  following  equation,  with  small  error. 
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0.943605ao  ♦  {ffifjy  -  -  «•***  » 

or  0.948605a0+  ysf'-f)  -  -  0.533667, 

or  0.948605a0  *  -  0.723907, 


whence  a  first  approximation  to  aQ  is 
a0(l)  =  -  0.76313 

The  superscript  -within  parentheses  is  used  to  designate  the  ordinal 
number  of  the  approximation. 

The  systematic  solution  of  (12)  continues  with  a  similar  modification 
of  the  first  two  equations  of  that  system.  Thus,  we  replace  the  first 
two  equations  "by  the  pair 


0.948605a0  -  0.944856a1  +  X  -X  I  „  -  0.583667, 

0.315529a0  ♦  0.944485a1  ♦  £  {  i  +  £J  -f  }  -  -  0.272021, 

1  m,  (i7) 

or  0.948605aQ  -  0.944856a1  =  -  0.614992, 

0.3l5529ao  +  o.944485ai  =  -  0.342141, 


which  yield  the  approximations 

aQ(2)  =  -  0.75718 

=  -  0,10930 

Application  of  this  technique  to  the  first  three  equations  of  (12) 


yields  the  system 


0.9^8605*0  -  0.9^85^  -  °-31«03a£  +  -  . 

0.315529^,  +  O.gyAS^  -  0.9y.2g^ 

O.I8915&V,  +  0.31*775^  +  O.^ltO^  *  f 
or, 

0.9^»8605a0  -  0.944856^  -  0.314805a,,  *  -  0.598655 
0.315529*0  +  0.9WJ85*!  -  0.9^225*2  «  -  0.293129 
0.18915680  +  0.314775*2  +  o.^i^oag  -  -  0.221696 

The  roots  of  this  system  are 

=  -  0.75352 
=  -  0.10684 
=  -  0.04823 

The  next  step  in  the  iterative  process  involves  the  solution 
system: 

0.948605ao  -  0.944856ax  -  0.314805*2  -  0.l8385la3  *  -  0.592820 

0.31552980  +  0.944485S2  -  0.944225*2  -  0.3]4710a3  =  -  0.283405 

0.18915600  +  0.314775*2  +  0.944l40a£  -  0.944073*3  =  -  0.192522 

0.13505aQ  +  0.l88847a2  +  0.3l4698a2  +  0.944o4la3  =  -  O.I63174 


(3) 
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(2) 


(1) 


0.583667 

-  0.272021 

-  0.176764 
(18) 


of  the 

(19) 


Solution  of  this  eysien  yields  the  approximations 


w 

V  =  -  0.75127 
<^(3)  «  _  0.105^9 
a^2)  *  -  0.04687 
a^1)  *  -  0.02865 


Hie  cert  step  in  ohe  refinement  of  the  approximations  to  the  a^'s 
entails  the  solution  of  the  following  system,  where  an  ellipsis  in  each 
equation  denotes  the  ommissioa,  for  typographical  reasons,  of  the  co¬ 
efficients  of  a^  and  &2  C  unchanged  from  system  (l9)J. 


0.948605aQ . 0,l88851a3  -  0.134882ai+  = 

0.31552980  +  .  .  .  -  0.3l4710a3  -  0.l888l6a4  = 

0.18915680  +  .  .  .  -  0.944073a3  -  0.31467984  « 

0.13505000  +  .  .  .  +  0„944o4la3  -  0.9440148^  = 

0.105009ao  +  .  .  .  +  0.3l4674a3  +  0.944001a4  » 


-  0.589984 

-  0.279434 

-  0.185904 

-  0.143320 

-  0.128733 


From  this  set  of  equations  we  deduce  the  approximations : 


a  <3)  „ 
2 

(2) 

a3  " 
a^^1)  = 


-  0.74978 

-  0.10463 

-  0,04610 

-  0.02779 

-  0.01953 


(20) 
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Tnrther  extension  of  the  procedure  leads  to  the  following  sys~em 
of  six  simultaneous  equations,  -which  ha3  "been  abbreviated  in  o.  manner 


similar  to  system  (20)'j 


0.948605a0 . 0.1B4882a^  -  0.104902a5  -  -  -  0. 588359 

0.315529a0  +  .  .  .  -  0.l888l6a^  -  0.134864a5  -  -  O.277345 
0.l89156a0  +  .  .  .  -  0.3l4679a^  -  0.l88802a5  -  -  0.182980 
0.13505080  +  .  .  .  -  0.9440148]^  -  0.3l4664a5  «  -  0.138447 
0.105009ao  +  .  .  .  +  0.944001aj!  -  0.94,398385  *  -  0.114113 
o.085900aQ  +  .  .  .  +  0.3l4664a^  +  0.943977^-  -  -  0.106107 


The  corresponding  values  of  the  a^s  are: 


(6) 

a0 

-  0.74871 

a<5>  = 

1 

-  o.io4o4 

- 

-  0.045590 

a3(3)  " 

-  0.027276 

a(2)  . 
a4 

-  0.018947 

.  U)  . 

-  0.014426 

Desk-machine  calculations  were  terminated  with  the  solution  of  the 
following  system  of  seven  simultaneous  equations: 


0.948605*0  -  .  .  . 

-  0.104902a5  -  0.085827*5  - 

-  0.587328 

0.315529*0  ♦  *  •  . 

-  0.134864a5  -  0.ic4892a^  - 

-  0.276085 

0.189156SQ  +  .  .  . 

-  0.lSS802a5  -  0.134857*5  - 

-  0.181360 

0.135050ao  +  .  .  . 

-  0.3l4664a5  -  0.18879735  « 

-  0.136179 

0.105009ao  +  .  .  . 

-  0.943983*5  -  0.314659*5  * 

-  O.IIO333 

O.0859OOSQ  +  .  .  . 

+  3.943977*5  -  0.943972ag  « 

-  0.094767 

0.072674aQ  +  .  .  . 

+  0.314656a,-  +  0.943967*^  » 

-  0.090127 

The  corresponding  roots  are: 


0^ 

n 

-  0.74793 

11 

-  0.10361 

*2(5)  - 

-  0.045235 

aW  . 

-  0.026940 

a^(3)  = 

-  0.018580 

a<2)  - 

a5 

-  0.013992 

U) 

a6 

~  0.011220 

Further  calculations  of  this  type,  involving  increasingly  larger 
systems  of  equations,  can  he  expeditiously  performed  hy  high-speed 
electronic  computers.  However,  the  preliminary  calculation  of  the 
numerical  coefficients  of  the  equations  is  more  formidable,  and  must 
necessarily  he  programmed  and  carried  out  as  a  separate  operation. 


(22) 
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The  relatively  elov  rate  of  convergence  way  be  deduced  frcn  the 
following  condensed  table  of  successive  fippraxiaaticcs  to  the  syetex  (12). 


TABUS  1 


n 

.(*) 

®o 

a1(n) 

a^Cn) 

»3(n) 

*5(a) 

a6(n) 

1 

-  0.7631 

-  0.1093 

-  0.0482 

-  0.0287 

-  0.0195 

-  0.0144 

-  0.0112 

2 

-  0.7572 

-  0.1C68 

-  0.0469 

-  0.0278 

-  0.0109 

-  0.0140 

3 

-  0.7535 

-  0.1055 

-  0.o46l 

-  0.0273 

.  0.0186 

4 

-  0.7513 

-  0.1046 

-  0.0456 

-  0.0269 

5 

-  0.7498 

-  o.io4o 

-  0.0452 

6 

-  0.7487 

-  0.1036 

7  -  0.7VT9 

The  convergence  of  each  sequence  [  \^n^}  i*  accelerated  by  applying 
the  Aitken  S^~  process  (6),  that  is,  the  associated  sequence  N(n)} 
is  formed  according  to  the  foraula 

i  (n)  .  4n~l)‘Sa(n*l)  -  «a(p)  (g, 

4“'1)  -  2ln(a)  + 

The  sequences  £  b^C11) ^  ,  m  *  0,  1,  2f  •  •  yield  the  following 
accurate  approximations  to  the  roots  of  the  infinite  system  (12): 


12 


»q  »  '  0-7^53 
a,  =  -  0.1C24 
ag  =  -  0.0444 
=  -  0.02c2 

84  =  -0.0179 
a_  =  -  0.0133 
ao  *  “  0.0105 


Ve  define  6  by  the  relation 

>W 

e  =  a  —  a. 

■Jyj,  “Vn.  >*V 

Id  the  particular  numerical  case  under  consideration,  €  assumes 
the  following  values,  deduced  frca  equation  (6)  in  conjunction  with  the 

preceding  set  of  final  approximations  to  the  a  *s. 

m 


6  -  -0.0535 
© 

et  =  +  0.0130 

6  =  0.0075 

e  =  0.0047 

£ 

€  -  0.0031 

-  0.0022 

^  =  0.0015 


The  £'s  thus  derived  are  then  used  to  evaluate- the  Kapteyn  series 
appearing  in  the  right  member  of  equation  (l).  To  this  end  ve  vrite 
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y*M> 


The  motivation  for  this  representation  of  the  Kapteyn  series  lies  In 
the  fact  that  the  series  comprising  the  first  term  of  the  right  side 
of  equation  (25)  has  been  evaluated  by  the  UEIVAD  for  p  «  3(1)6, 

/JL  »  0.25(0.25)2.0(0.5)6.0,  and  yfa  -  0.2(0.1)0.8(0.05)0.95,  O.975. 


The  multiplicative  factor  appears  In  equation  (25)  because 

the  ultimate  goal  of  the  calculation  Is  the  evaluation  of  the  Goldstein 
factor  K.  ,  vhere,  as  already  stated  In  equation  (3), 


H-j*?'  fV u) 
ZXAaTIT 


The  effect  on  K  of  introducing  the  correct  coefficients  a^  In 
place  of  the  approximate  coefficients  affl*  may  be  inferred  from  an 
examination  of  the  following  table. 


TABLE  2 


Qoldstein  Factor  K  for  p  *  3  and  ■  1.5 


Ms, 

(ayprox.  value 
based  on  a^*) 

/^(accurate  value 
based  on  as) 

0.2 

1.4268 

1.4034 

0.3 

1.0972 

1.0756 

0.4 

0.9163 

0.8945 

0.5 

0.7964 

0.7731 

0.6 

0.7016 

0.6760 

0.7 

0.6123 

0.5837 

0.0 

0.5124 

0.4808 

0.85 

0.4516 

0.4191 

0.9 

0.3766 

0.3442 

0.95 

0.2723 

0.2428 

0.975 

0.1932 

0.1678 

The  relative  error  in  the  above  values  of  %  based  on  the 
approximate  coefficients  a^*  increases  from  about  2  percent  to  more 
than  15  percent  as  approaches  unity.  On  the  other  hand,  the 

absolute  error  in  the  approximate  values  of  remains  fairly  constant, 
varying  from  2.16  X  10“^  to  3«25  X  10”^. 
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Appended  to  this  report  are  tables  of 


•f 


i  +  U  i 


>n-«  "  (Xr^+l)*- 


and 

nf  A  Z  < 

^  r  V^+4)"f  ^ 


derived  from  similar  results  computed  by  UNIVAC  in  the  solution  of 
Applied  Mathematics  Laboratory  Problem  42-54. 

These  data  can  be  used  as  the  basis  for  further  computation  of 
accurate  values  of  K  in  the  manner  illustrated  in  this  report. 

An  investigation  was  also  made  of  the  applicability  of  published 
theorems  gf  von  Koch^)  and  March  (8) 

to  the  problem  of  determining  the 

existence  of  a  solution  to  the  infinite  system  (2).  The  sufficient 
conditions  contained  in  these  theorems  are  too  restrictive  to  be  satisfied 
by  the  system  under  consideration  in  this  report,  and  consequently,,  the 
theorems  do  not  apply.  A  rigorous  proof  of  the  existence  of  a  solution 
of  a  system  such  as  (2)  apparently  is  not  yet  available. 

I  should  like  to  conclude  vith  an  acknowledgment  of  the  material 
assistance  rendered  by  Mrs.  Anita  B.  kllam  by  her  extensive  calculations 
of  the  modified  Bessel  functions  of  half-odd-integer  orders,  which 
entered  into  the  evaluation  of  the  coefficients  of  the  numerical  system 
of  equations  treated  in  this  report. 
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